Mosaic multi-state scenario vs. one-state description of supercooled liquids 
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According to the mosaic scenario, relaxation in supercooled liquids is ruled by two competing 
mechanisms: surface tension, opposing the creation of local excitations, and entropy, providing the 
drive to the configurational rearrangement of a given region. We test this scenario through numerical 
simulations well below the Mode Coupling temperature. For an equilibrated configuration, we 
freeze all the particles outside a sphere and study the thermodynamics of this sphere. The frozen 
environment acts as a pinning field. Measuring the overlap between the unpinned and pinned 
equilibrium configurations of the sphere, we can see whether it has switched to a different state. We 
do not find any clear evidence of the mosaic scenario. Rather, our results seem compatible with the 
existence of a single (liquid) state. However, we find evidence of a growing static correlation length, 
apparently unrelated to the mosaic one. 
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It is common opinion that in finite dimension a diver- 
gence of a relaxation time r at nonzero temperature is 
associated to a diverging characteristic length £. The 
idea is that when this length increases, relaxation pro- 
ceeds through the rearrangement of ever larger regions, 
taking a longer and longer time. The relation between r 
and £ depends on the physical mechanism of relaxation. 
Two main mechanisms are activated relaxation of a ip- 
dimensional droplet of size £, giving r ~ exp(A£^/T), 
and critical slowing down, where r ~ £ z []]]. 

Glass-forming liquids are tricky: relaxation times grow 
spectacularly (more than ten decades) upon lowering the 
temperature, without clear evidence of a growing static 
cooperative lenght. In particular, density fluctuations are 
thought to remain correlated over short distances close 
to the glass transition (although there are some indica- 
tion that energy fluctuations might develop larger cor- 
relations @, Q). Thus the concept of dynamic hetero- 
geneities is central to several theories of the glass tran- 
sition @, [E IE 01' where the role of order parameter is 
played by dynamic quantities such as local time corre- 
lators, which become correlated over the growing dy- 
namic lenght scale £dyn- No thermodynamic singularity is 
present in these theories. Dynamic singularities are also 
typically absent at finite temperatures, with the notable 
exception of mode coupling theory (MCT) Q, recently 
recast in terms of dynamic heterogeneities M. Note, how- 
ever, that the experimental values of £dyn fiol [ill [l3j 
are barely in the nm range, the same as density correla- 
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The mosaic scenario (MS) [15 . 
the conceptual framework of nucleation theory, identifies 
on the other hand a static correlation length. Deeply 
rooted in the physics of mean-field spin glasses, the MS 



crucially assumes the existence of exponentially many 
inequivalent states exp(iVE), below the mode coupling 
temperature Tmc (S is called complexity or configura- 
tional entropy, and N is the size of the system). Suppose 
the system is in a state a and ask: what is the free en- 
ergy cost for a region of linear size R to rearrange into 
a different state /3 with the same free energy? Accord- 
ing to the MS, there are two opposing contributions: a 
surface cost due to the mismatch of a and /3, propor- 
tional to TR e , where T is a generalized surface tension 
(d < d— 1), and an entropic gain proportional to TY,R d , 
arising from the fact that the larger the region, the higher 
the number of possible rearrangements. This is similar 
to nucleation theory, with TS playing the role of the 
free energy difference between the two phases. For small 
droplets, the surface contribution dominates, whereas for 
large R the volume contribution wins, and eventually the 
rearrangement occurs. Thus, similarly to nucleation the- 
ory, a length scale emerges, £ mos = (T/TS) 3 ^ fixed by 
the balance of the two contributions. A droplet with 
R < 6nos may flip to a different state, but the pinning 
field provided by the surrounding system will make it 
flip back; for R > £ mos the region is thermodynamically 
favoured to flip, but the larger R the longer the time 
the rearrangement takes. £ mos thus acts as the typical 
cooperative length. Furthermore, by noting that £ mos 
diverges at the temperature T$ where £ vanishes, the 
MS provides a rationale for the Vogel-Fulchncr-Tamman 
(VFT) law for the relaxation time, r ~ exp[A/(T - lb)]. 

In this Letter we report a numerical test of the MS. 
For this purpose, the formulation of ref.[l8|is particularly 
convenient. Imagine picking a reference equilibrium con- 
figuration and freezing all particles except those within a 
sphere of radius R (containing M particles). This region 
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is thus embedded in a very large box of frozen particles 
which act as a pinning field. The key point of the MS 
is that a sufficiently large sphere will be thermodynami- 
cally favored to flip to a different state. If the reference 
configuration is in state a, and calling (3 one of the ex- 
ponentially many different states, the MS gives [l8| 
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P aa (R) = £ expfW/T), Pa0 (R) = - cxp(£i? d ), 

(1) 

as the probabilities for the sphere to remain in the ref- 
erence state a and to flip to f3, respectively (the normal- 
ization is Z = exp(Ti?7T) + exp(£i? d )). The mosaic 
fragmentation of a state into regions of size £ m os corre- 
sponds to an exponentially sharp jump from p aa ~ 1 to 
Paa ~ at R ~ £ mos [18J . Let us call the configura- 
tion of the sphere thermalized with a given environment 
the pinned equilibrium state of the sphere. If we define 
a suitable overlap q to measure the similarity in phase 
space between this pinned state and the reference state 
a, the MS predicts, 
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where q aa is the self-overlap of state a and go is the 
typical overlap between different states. Thus within the 
MS one expects a first order transition [191 ] , where qms (R) 
drops (or at least makes a crossover) at £ mos from q aa 
to go- Note finally that for large R the overlap decays 
exponentially, qus(R) ~ Qo ~ exp[-(£i? d - YR e /T)\. 

We have realized this gedanken experiment for the soft- 
sphere binary mixture [20(, a simple fragile glassformer 
(see ref. [2l| for an attempt to find a mosaic length-scale 
in spin models using this formulation). We used the swap 
Monte Carlo algorithm [22] (for simulation details see 
ref. 1231). Equilibration time is considerably shortened, so 
we can perform equilibrium simulations of large systems 
well below T MC =0.226 [HHl- This is important, since 
in the standard MS Tmc acts like a spinodal tempera- 
ture 15j , above which there are not many states and the 



surface tension is zero [3l| . We first equilibrated systems 
both of 2048 and of 16384 particles, then picked reference 
configurations used to equilibrate M particles within of 
the sphere while keeping the remaining frozen. We stud- 
ied M = 5, 20, 50, 100, 200, 400, 800, 1600, 3200, and 
5500 at temperatures T/T MC = 4.42, 2.13, 1.54, 1.15, 
0.94, and 0.89. The largest R was 10.95, one order of 
magnitude larger than the particle size. Though remi- 
niscent of the work of ref. l25l . we focus on static rather 
than dynamic quantities. 

To introduce the overlap between the reference state a 
and the pinned state of a sphere of radius R, we divide 
the space in cells of side I and define, 
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where rij is the occupation number of cell i. The sum runs 



FIG. 1: Overlap q(R) for T/Tmct =4.42 (circles), 1.54 
(squares), 0.94 (triangles), 0.89 (inverted triangles). The 
largest value of R corresponds to M = 5500 particles in the 
sphere. Lines are the fit from the IS argument (Eq.[6]). Inset: 
data for all temperatures scaled according to Eq. (0. 



over all the cells in the sphere and M = 4/3irR 3 p. Oc- 
cupation numbers are averaged over many independent 
configurations, both of state a (between 4 and 16) and 
the pinned state (10 to 100). The overlap of two identi- 
cal configurations of the sphere is q — 1 , whereas for two 
independent configurations q = qo = e, with e = pi 3 . I is 
such that e <C 1, but larger than the typical vibrational 
amplitude of the particles [23|. Here e = 0.06. 

Our results for q(R) are shown in Fig. [T] at four differ- 
ent temperatures, two above and two below Tmc- The 
first feature we notice is that q(R) is always a smooth 
function, with no clear jump, nor evident crossover val- 
ues of R. Yet, even at the lowest temperature, q becomes 
as low as 0.2 at the largest value of R, showing that the 
sphere has in fact largely decorrelated with respect to the 
reference state a. This means that whatever the correct 
relaxation mechanism at low temperatures, this mecha- 
nism is at work here, since on average 80% of the parti- 
cles have rearranged at a temperature well below Tmct- 
Moreover, q(R) does not seem to decay exponentially to e 
for large R, but rather as a power law (see later). Finally, 
there is no clear indication of a plateu corresponding to 
the self-overlap q aa . Unsurprisingly, a mosaic fit of the 
data, using Eqs. [1] and [2l gives unphysical values of the 
paramenters, in particular a negative value for the com- 
plexity S and of the surface tension T. 

Three objections can be raised at this point. First, 
the spheres are not large enough, the mosaic drop takes 
place for £ mos larger than our largest R. To this there 
are two replies: first, £ mos should decrease for increasing 
T, so it is unclear why we do not see anything while ap- 
proaching Tmc from below; second: as already said, the 
largest sphere has largely decorrelated, so, if the mosaic 
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bit of decorrelation still has to come, it will contribute 
only to 20% of the rearrangement. Second objection: the 
equilibration time of the sphere is longer than our runs, 
and our results for q(R) are obtained in a metastable re- 
gion. To check this we have repeated the semulations but 
initializing the sphere in a different equilibrium state (3, 
having zero overlap with the reference state a. What we 
find (Fig. right) is that the overlap as a function of 
time increases, thermalizing at the same value as when 
starting within the same state a. Thus, even for the 
largest R, the asymptotic value of q(R) does not depend 
on the initial configuration of the sphere, and there is no 
hysteresis. Third objection: the lowest temperature is 
not low enough, it is too close to Tmc, an d the mosaic 
mechanism is not yet at work. We cannot exclude this. 
However, at our lowest temperature standard molecular 
dynamics is completely stuck, so one would expect the 
MS to describe the relaxation. Besides, an estimate of To 
from a VFT fit of the relaxation time gives To « 0.80Tmc 
(see Fig. [21 left), so one would really expect the mosaic 
mechanism to be operating at these temperatures. 

If we give up the key assumption of the MS, that is 
the existence of many states at low temperature, we can 
interpret our data as relaxation within a single state with 
self-overlap equal to e, i.e. the liquid state. The argument 
goes as follows. Let us divide the sphere in shells of radius 
r. From Eq. [3] we have, 
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where G(r) = (q(r))/e is the average overlap per unit 
volume between the state a and the asymptotic pinned 
state, at distance r from the centre of the sphere. The 
effect of the pinning border at the centre of the sphere is 
expected to decay as exp(— i?/A), where A is a correlation 
length 26] . If we say that the unpinned state (the liquid) 
has self-overlap e, and we assume G{r = R) = 1 (sticky 
boundary), a reasonable form for G(r) is given by, 

G(r) = (l-e) e^ R - r ^ x + e. 

Plugging G(r) into Eq. [4] and having defined x 
we obtain the one-state (IS) overlap, 
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This function describes the relaxation within a single 
state with self-overlap e, and it is quite different from 
the MS form (Eq. [2|). In particular, qi S (R) - e ~ \/R, 
for large R. In Fig. [1] we report the fit to the data ob- 
tained with the IS overlap, which seems quite reasonable 
[32I ] . The only fitting parameter is the correlation length 
A, which increases by a factor 7 with decreasing T in our 
temperature span (Fig. [3]). 

We can test assumption (J5J) by studying, at fixed R, 
the average overlap as a function of the distance r from 
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FIG. 2: Left: Relaxation time vs. temperature for MC swap 
dynamics and VFT fit. r is an integrated relaxation time ob- 
tained from the coarse-grained density autocorrelation stud- 
ied in ref. [2^ . Right: Instantaneous overlap vs. time for 
a and f3 initial configurations of the sphere (T = 0.89Tmc, 
M = 5500). 



the center (Fig. [3J inset). We plot log[(G(r) - e)/(l - e)] 
vs. R — r, for two different values of R, at our lowest 
T. The data are compared to our guess Eq. [5l with A 
obtained from the fit of q(R). The agreement is reason- 
able. Another test is to plot [qis(R) — e]/(3 — 3e) vs. R/X, 
where according to © data for all temperatures should 
fall on a single master curve. Again (Fig.[l] inset) we find 
a nice agreement. We note that the correlation function 
G(r) is a kind of the "point-to-set" correlation functions 
described in [27]]. Other than testing the MS, it seems 
from our results that G(r) may be a useful tool to detect 
growing spatial correlations in supercooled liquids. 

The one-state decay of q(R) described above should in 
fact be used also in the small R regime of the mosaic 
scenario. This modifies Eq. [3 as follows: 



'ims (R) (R) qis(R) +p a p(R)t 
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where now the limiting self-overlap of q±s(R) for R — > 00 
is q aa , hopefully larger than e in a multi-state scenario. 
Thus, the correct mosaic prediction is that the relaxation 
of the sphere within state a is ruled for small R by qis (R), 
but it is interrupted at R = £ mos by the entropic drive 
kicking in, such that q(R) — > e for R > £ ms . We have 
tried Eq.[7|to fit the data of Fig. [2 leaving 8 fixed to avoid 
having too many parameters. For 8 = 2 and 8 = 3/2 
(two values one might expect [Hj]), the fits locate the 
mosaic drop just beyond the largest simulated size, with 
qaa ~ £■ This makes the mosaic fit rather awkward. In 
our opinion, given the present data, it is more natural 
to conclude that a is the (unique) liquid state, with self- 
overlap e, and that there is no mosaic drop beyond our 
largest R. 
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FIG. 3: Correlation length A vs. T. Inset: [G(r) - e]/(l - e) 
vs. R-r for T = 0.89T M c at M = 400 and M = 5500. The 
line is our prediction, Eq. [5] 



We have studied the themodynamics of a sphere em- 
bedded in a frozen equilibrium environment. Our data 
are compatible with the naive expectation that the in- 
ner shells of the sphere decorrelate and rearrange more 
than the outer shells close to the border. The penetra- 
tion scale of the pinning border is given by a length A, 
which increases sharply with decreasing temperature. In 
the paramagnetic phase of the Ising model, A diverges 
for T — > T e , just as the standard correlation length [28j |. 
This means that A has physical significance and begs for 
an investigation of the relation between A and r also in 
glass-formers. Such an investigation, however, would re- 
quire a finer resolution of A(T) around and below XjyiCj 
which we leave for a future study. The length A is neatly 
defined through the novel correlation function G(r). Our 
study is purely thermodynamic, and for this reason we do 
not see any crossover at the mode coupling temperature. 

Our data do not show any clear evidence for the mosaic 
mechanism. There is the possibility that the experiment 
or the observable considered are not the most suitable. 
For example, if the mosaic excitations are highly non- 
compact or fractal [2^ | , it may be argued that a different 
experiment would be needed. Having said that, we re- 
main with the impression that the one-state scenario fits 
reasonably well our data. Either there is no surface ten- 
sion, or it is very weak and diluted over a wide region 
across the border. This suggests that, if a multi-state 
scenario is valid, a version more sophisticated than the 
here tested is needed. This novel description should not 
only reproduce the numerical data, but fit them signifi- 
cantly better than the single-state scenario. 
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